Magnetization plateau and incommensurate spin modulation in Ca3Co206 
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The magnetic properties of a trigonal prism unit of the spin-2 frustrated compound Ca3Co206 are studied by 
means of the density-matrix renormalization group method. A magnetization plateau at ms/3 (rua is the satura- 
tion magnetization) with ferrimagnetic structure is observed. By fitting the experimental data of magnetic curve, 
an estimation of the couplings gives Ji=— 26.84K, J2=0.39K, and J3=0.52K. The local magnetic moments are 
unveiled to exhibit an incommensurate sinusoidally modulation along the three chains of the trigonal prism, 
which gives a strong theoretical support to the experimentally observed incommensurate partially disordered 
antiferromagnetic state for Ca3Co2 06. The present result suggests that the modulation indeed originates from 
the competition of antiferromagnetic and ferromagnetic couplings. 



Frustrated quantum magnets have been extensively studied 
owing to their exotic properties |2]. Among others, the first 
experimental realization of feiTomagnetic (FM) chains cou- 
pled with antifeiTomagnetic (AFM) interactions in a triangu- 
lar plane in the spin frustrated chain compound Ca3Co206 has 
attracted much attention in the last decade |3,-^. Ca3Co206 
consists of spin chains with alternating face-sharing octahe- 
dral (Co I) and trigonal prismatic (Co II) CoOg polyhedra 
along the c axis, and aiTanged in a triangular lattice in the ab 
plane |3], as shown schematically in Fig. 1(a). The Co'^+ ions 
on the Co I and Co II sites have low-spin {S=0) and high-spin 
(S=2) states, respectively|6]. Within the chains, direct Co-Co 
overlap leads to a strong FM intrachain interaction Ji (about 
25K) 1 7]. For the AFM interchain super-superexchange inter- 
actions J2 and J3, as the overlap of the O 2p orbital is very 
small, J2 is weaker than J3 |l7|]. 

A number of techniques (x-ray ^ |^, neutron scattering 
ifToll . NMR lH], calorimetry iH, and fiSR fH, etc.) have 
been used to study this intriguing compound in different as- 
pects, such as the nature of the order at different temperatures, 
the magnetization process under various conditions, and so 
on. It was found that the compound has two transition tem- 
peratures Tcic±25K and Tc2—7K, between which the magne- 
tization curve has a TOs/3 plateau (m^ is the saturation mag- 
netization) [3]. Below Tc2, the magnetization shows a mul- 
tistep structure and the hysteresis, which was suggested in a 
spin freezing state IHQl, and appears still under debate. An- 
other intriguing issue is the magnetic state between Td and 
Tc2 in absence of a magnetic field. Owing to the spin frus- 
tration, two possible magnetic structures were proposed: (1) 
a ferrimagnetic (FI) structure, with two third of spins up and 
one third spins down; (2) a partially disordered antifeiTomag- 
netic (PDA) structure, with one third spins up, one third spins 
down, and one third spins incoherent. 

Recently, the zero-field magnetic state of Ca3Co206 be- 
tween Tci and Tc2 has been shown to be a sinusoidally modu- 
lated incommensurate PDA (iPDA) state with a modulation of 
very long periodicity (^lOOOA) along the chain direction by 
resonant X-ray scattering [Ist] and neutron diffraction iH 
measurements. Such a magnetic structure cannot be explained 
in a simple two-dimensional (2D) triangular Ising model JTsI, 
in contrast to some previous works. Therefore, to understand 
the iPDA state, the trigonal prism unit of Ca3Co206 includ- 



ing allpossible Heisenberg interactions should be considered 
lll4l[l6ll . In the trigonal prism, the AFM interchain couplings 
J2 and J3 follow helical paths and connect the Co II (5=2) 
sites in adjacent chains [Fig. \lla)]. The competition between 
the AFM and FM interactions was argued to be responsible 
for the modulation [14. ,16.1 . The emergence of the iPDA state 
was suggested to be related to the trigonal prism unit lfl6ll . 
However, a direct theoretical analysis of the iPDA state in the 
trigonal prism unit of Ca3Co206 is still absent. 
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FIG. 1: (Color online) (a) Schematic structure of the trigonal prism 
unit of Ca3Co2 06, Ji (black lines) FM; J 2 (blue dotted lines) AFM; 
J3 (red dashed lines) AFM. Co II (S=2) sites in different chains 
are colored differently, (b) Sketch of the equivalent S=2 spin chain 
model [Eq. ([TJ]. 

In this paper, by means of the density-matrix renormaliza- 
tion group (DMRG) method we shall study the magnetic prop- 
erties of the 5=2 spin chain model with the trigonal prism unit 
under various circumstances, and fit the experimental result of 
magnetization curve of Ca3Co206 to estimate the three dif- 
ferent couplings Ji, J2 and J3. The origin of the iPDA state 
in Ca3Co206 will also be addressed. Our results show that 
the system exhibits a ms/3 magnetization plateau with a FI 
structure, and a very good fitting to the experimental data is 
obtained, with the theoretical fitting values in well agreement 
with those obtained by other experimental means. The in- 
commensurate sinusoidally modulated structure with a long 
period along the chain direction is observed in the local mag- 
netic moments in the present model, which may shed light on 
the origin of the iPDA state in the compound Ca3Co206. 

The trigonal prism unit of Ca3Co206 [Fig. 1(a)] can be 
equivalently expanded to a S=2 Heisenberg spin chain with 
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FIG. 2: (Color online) Magnetization curves of spin-2 trigonal prism 
chain with (a) J2/I Ji|=0.2 for different J3, and (b) J3/I Ji|=0.2 for 
different J2. Magnetic phase diagrams in the h — J2,z plane of the 
S=2 trigonal prism chain under interest for (c) J2 / 1 Ji | =0 . 2, and (d) 

J3/|Ji|=0.2. 



nearest neighbor, next, next next, and next next next nearest 
neighbor interactions, as illustrated in Fig. [lib). The Hamil- 
tonian under interest is given by 

L/3 

H — ^^[A(S3i-2 • 832+1 + Sa.i-i • 83^+2 + 83^ • 83.^+3) 

i=l 

+J2{^d,i-1 • 83^+1 + 83; • 83^+1 + 83^ ■ 83i-|_2) 
+■^3(83^^2 • 83.i_i + 83.i_i • 83; + 83i_3 • 83^+1)] 
L/3 

- hY^Sh (1) 

where 8j is the S=2 spin operator at the y-th site, L is the to- 
tal number of spins in the trigonal prism, yi<0 (/2,3>0) is 
the FM (AFM) interaction, and h is the magnetic field applied 
along the chain direction. For convenience, we take q^B=^- 
During the DMRG \ Yf\ iteration, the length of the chain is 
taken at least 90, the optimal states are kept 120, and the trun- 
cation error is less than 10^^. 

Figs- Eta) and (b) show the DMRG results of the mag- 
netization curve of the 5=2 trigonal prism spin chain sys- 
tem for different couplings. It can be seen that under proper 
couplings, the system has a m^/S magnetization plateau 
that is consistent with the necessary condition |18] n{S — 
m)=integer for the appearance of magnetization plateau in 
Heisenberg quantum antiferromagnets, where n is the period 
of ground state, S is the magnitude of spin, and m is the 
magnetization per site. For the present system, n=3. In the 
plateau state, the local magnetic moments are found to have 
the same value in each chain of the trigonal prism, arranged 
in a FI structure, namelyf— 2, 2, 2), in agreement with the ex- 
perimental observation yfl. The width of the plateau increases 
with increasing AFM interactions J2 and J3, and decreases 
with increasing FM coupling | Ji|, which are consistent with 
the behaviors of magnetization plateau observed in other sys- 
tems 1 19]. Through massive calculations, the magnetic phase 
diagrams in the h — J2( J3) plane can be obtained by observ- 
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FIG. 3: (Color online) The experimental data of magnetization curve 
for CasCozOs at r=10K from Ref. [H are fitted with the DMRG 
results, giving Ji=-26.84K, J2=0.39K and J3=0.52K. In the cal- 
culation, the number of Co sites is chosen as 180. 



ing the critical and saturation magnetic fields of magnetization 
curves under different couplings, as depicted in Figs. [Jfc) and 
(d). It appears that the region for the appearance of ms/3 
magnetization plateau is narrow. 

The nis / 3 magnetization plateau of the trigonal prism chain 
compound Ca3Co206 has been already observed experimen- 
tally llT2ll . To fit the experimental data, we make a lot of 
DMRG calculations ll2lll. The fitting results for the data 
of T=10K from Ref. are presented in Fig. [3] which 
shows a good agreement, giving Ji=— 26.84K, J2=0.39K, and 
J3=0.52K. The so-obtained Ji is very close to 25K estimated 
from the experimental study ^ and first-principle calcula- 
tions |7] of Ca3Co206, and Js/Ji^— 0.019 is also consistent 
with the value —0.018 obtained in Ref. within a classical 
treatment. Below the plateau, as shown in Fig. [3] the mag- 
netization curve increases very sharply as indicated by a cusp 
observed experimentally |20] in the magnetic susceptibility, 
which was once regarded as a feature of FI state of Ca3Co206 
between Td and Tc2- Our DMRG results indicate that the 
iPDA state can also give rise to a similar behavior below the 
magnetization plateau. 

Now let us look at the incommensurate magnetic modula- 
tion observed in Ca3Co206 through the trigonal prism spin 
chain model [Eq. ([TJ]. For this purpose, the magnetic mo- 
ments are calculated. Figure |4] presents the DMRG result 
of local magnetic moment (Sj) for Ji:J2:J3=— 1:0.3:0.4 as 
an example. To see clearly the magnetic modulation, {S^j), 
(•^Ij-i) ^^'^ {^3j-2)' which coiTespond to the three chains of 
the trigonal prism in Fig. [Tfa), are studied, respectively. It is 
evident that the spatial dependence of (Sj) for the three chains 
along the c axis exhibits modulations sinusoidally around 
0.986 ({Slj)), 0.414 and -1.4 ((5|,_2)), indicat- 

ing that the state is neither FI nor PDA, but an iPDA, thereby 
supporting the recent experimental observations |9, 13, 14]. 
During the DMRG calculations, the systems with various 
chain lengths are compared. It is found that the expectation 
value of each local magnetic moment (S^) is almost indepen- 
dent of the system size, indicating that the modulation is not 
a finite-size effect. In Fig. ID we show the case with 120 
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FIG. 4; (Color online) The spatial dependence of the local mag- 
netic moments (Sj) of the S=2 trigonal prism spin chain model for 
Ji:J2:J3=-l:0.3:0.4. 



sites (40 unit cells) of the S=2 spin chain [Eq. (1)] that gives 
the period of modulation about 36 unit cells, showing that the 
system can have a magnetic modulation with a long period. 
Besides, we have also uncovered the less |Ji| (|Ji|>J2,3), 
the shorter the modulation period. When J2 or J3 is set zero 
and I Ji I dominates, the similar magnetic modulations are also 
seen. In a general case with | Ji |> J2,3, the period of modula- 
tion is observed to decrease with increasing AFM interactions 
J2^3 and increase with increasing FM coupling | Ji | . These ev- 
idences show that the modulation indeed originates from the 
competition between the FM and AFM interactions. 

It was argued that the magnetic modulation in Ca3Co206 
may be caused by both the single-ion anisotropy and the com- 
petition of couplings 1 14]. To clarify this point, we consider 
the effect of single-ion anisotropy —Dj^iiSi)^ iD>0) on 



the magnetic properties of the spin-2 Heisenberg complex 
chain system [Eq. (1)]. The results show that the single-ion 
anisotropy weakens the modulation and prolongs the modu- 
lation period. For instance, when D=0.2, the period exceeds 
50 unit cells. With further increasing D, the modulation van- 
ishes and the local magnetic moments in different chains form 
a FI structure, i.e., (—2, 2, 2) when DiOA. When we replace 
Ji terms in Eq. (1) along the three chains as Ising couplings 
( Ji> J2,3), we find no modulation for the magnetic moments, 
implying that the 2D triangular Ising treatment cannot de- 
scribe the observed iPDA state below the plateau. 

In summary, the magnetic properties of the 5=2 trigonal 
prism unit of Ca3Co206 are studied by means of the DMRG 
method. The ms/3 magnetization plateau with a FI struc- 
ture is observed, which is consistent with the experimen- 
tal observation. The phase diagrams in the h — J2 3 plane 
are obtained. The experimental data of magnetic curve of 
Ca3Co206 are nicely fitted to the DMRG results, giving an 
estimation of the interactions Ji=— 26.84K, J2=0.39K, and 
J3=0.52K, in agreement with the previous experimental mea- 
surement. The iPDA state with a long period observed experi- 
mentally is supported by the calculated incommensurate sinu- 
soidally modulated magnetic moments in the trigonal prism 
unit of Ca3Co206. The effects of single-ion anisotropy and 
Ising limit are also discussed. Our calculation gives a strong 
support to the proposal that the iPDA state originates from the 
frustration due to the competition between the intrachain FM 
interaction and the interchain AFM couplings, and indicates 
the crucial role of the trigonal prism structure in understand- 
ing the iPDA state of Ca3Co206. 
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